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Abstract. We present results from numerical simulations of nonlinear MHD dynamo action produced by three- 
dimensional flows that become turbulent for high values of the fluid Reynolds number. The magnitude of the 
forcing function driving the flow is allowed to evolve with time in such way as to maintain an approximately 
constant velocity amplitude (and average kinetic energy) when the flow becomes hydrodynamically unstable. It 
is found that the saturation level of the dynamo increases with the fluid Reynolds number (at constant magnetic 
Prandtl number), and that the average growth rate approaches an asymptotic value for high fluid Reynolds 
number. The generation and destruction of magnetic field is examined during the laminar and turbulent phase of 
the flow and it is found that in the neighborhood of strong magnetic flux "cigars" Joule dissipation is balanced 
by the work done against the Lorentz force, while the steady increase of magnetic energy occurs mainly through 
work done in the weak part of the magnetic field. 
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1. Introduction 

The dynamo problem is mainly related to the understand- 
ing of the processes of magnetic field generation due to the 
motions of conducting fluids, the so-called dynamo action 
(CowlingUHMI Moffatt 119751 Parker [1375) 1 . One often di- 
vides dynamo theory into two theoretical regimes: First, 
the kinematic dynamo problem in which magnetic fields 
and flows are decoupled. The velocity of the flow is pre- 
scribed and the evolution of the magnetic field is governed 
by the induction equation (here in dimensionless form): 
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where B is the magnetic field and u is the prescribed ve- 
locity field. The quantity Re m = u£/r) is the magnetic 
Reynolds number and it is often huge in astrophysical sys- 
tems, where u and I are characteristic velocity and length 
scales and rj is the magnetic diffusivity. A flow acts as a 
fast dynamo if the growth rate does not tend to zero in 
the limit of infinite magnetic Reynolds number. The so- 
lar dynamo is an astrophysical example of fast dynamo 
action since it operates on the convective time scale of 
the fluid which is fast compared to that of the Ohmic dif- 
fusion. In the fast kinematic dynamo problem one seeks 
velocity fields that lead to the exponential amplification of 
the magnetic field through induction Eq. lfT)l. A powerful 
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mechanism for growth of the magnetic field is the stretch- 
ing of the magnetic field lines (cf. Archontis ct al. 2003), 
and is primarily achieved by chaotic flows. The kinematic 
approach is valid when the magnetic field is weak and 
there are no further dynamical effects. In the second, fully 
dynamical dynamo regime, the velocity field is not pre- 
scribed and the exponential growth of the magnetic field 
saturates when the field becomes strong enough to modify 
the flow sufficiently through the feedback by the Lorentz 
force. The saturation state of the dynamo then emerges 
as a self consistent solution to the equations of nonlinear 
magneto-hydrodynamics. In a fully self-consistent dynam- 
ical experiment the forcing would have internal causes, e.g. 
of convective origin — here we are concerned with the case 
of an externally prescribed forcing function. 

The nonlinear properties of fast dynamos have received 
considerable attention and a variety of both analytical and 
numerical studies have provided a valuable insight into the 
nature of the dynamo action and saturation (Nordlund et 
al. IT5521 Brandenburg et al. IT5551 Cattaneo et al. IT55TH 
Zienicke et al. H998l Brummell et al. H998[l . However, a dif- 
ficulty that arises in numerical simulations is that often, 
fast dynamo flows become hydrodynamically unstable at 
fluid Reynolds number (Re= ui/v, v being the viscosity) 
greater than a critical value (Podvigina & Pouquet 1994) 
and spatio-temporal turbulence appears to be nascent for 
higher values of Re. As a result the desired flow velocity 
is modified even in the kinematic regime and the problem 
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contains the complexity of that of a turbulent dynamo. 
Thus, the choice of the forcing function that drives the 
flow is crucial. We use a forcing function with an ampli- 
tude that is allowed to evolve with time, keeping the aver- 
age kinetic energy approximately constant through both 
laminar and turbulent phases. The growth rates and the 
saturation level of such turbulent dynamos are examined 
with respect to the magnetic and fluid Reynolds num- 
bers. Also of great interest is the understanding of the 
processes at work in the kinematic (both laminar and tur- 
bulent phase of the flow) and in the saturated regime. We 
present an analysis of the Lorentz work and Joule dissi- 
pation, which shows that they are in close balance and 
that dynamo-work occurs mainly in the weak part of the 
magnetic field where dissipation is insignificant. 

The paper is organized as follows. Section [21 contains 
the equations, the numerical method and a description of 
the driving of the flow that has been used in the simula- 
tion. Results for turbulent dynamo action are presented 
in Section |3J The structures of the velocity field and the 
magnetic field in physical space are discussed in Section^] 
Whether or not the growth rates and the saturation level 
of the dynamo depend on the magnetic or fluid Reynolds 
number is examined in Section and Section El A quali- 
tative understanding of the nature of the dynamo is dis- 
cussed in Sectiond using results from a kinematic dynamo 
experiment. Section |S1 contains the overall conclusions of 
the present numerical simulations. 

2. The simulation 

2.1. The equations 

The compressible MHD equations are solved numerically 
in a periodic computational domain, with periodicity of 2ir 
in all three directions. Apart from the magnetic induction 
equation Eq. JTJ, the additional equations solved are the 
following: 



dp 

at 

dt 
de 

dt 



-V ■ pa, 

-VP + j x B + f - V • (puu), 
-V • (eu) - PV • u + Q v + Qj 



^cool 



(2) 
(3) 
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where p is the fluid density, P the pressure, j is the electric 
current density, e is the internal energy and f is an exter- 
nal forcing term. Q v and Qj are the viscous and Joule 
dissipation respectively. In the experiments, a Newtonian 
cooling term, Q coo i = (T— Xo)/ T cooi is used, where T is the 
temperature and t coo i = 5. The cooling term only works 
as a sink of the thermal energy, in order to balance the 
heat generation by dissipative terms in the energy equa- 
tion and the actual magnitude of the cooling time is of 
little importance with no effect on the results of our sim- 
ulation. 

The above equations are solved numerically on a stag- 
gered mesh. The time stepping is performed by a third 



order predictor-corrector method (Hyman I1979f) . The 
derivatives and interpolations are of 6th and 5th order re- 
spectively and the numerical scheme conserves V • B = 
exactly. Numerical solutions are obtained on a grid of a 
maximum of 180 3 points, using a modified version of the 
code by Galsgaard, Nordlund and others (Galsgaard & 
Nordlund HnS3 Nordlund et al. 11555)1 . The initial mag- 
netic field is chosen to be a weak random perturbation 
with an amplitude of 10~ 5 in non-dimensional units. 

2.2. Flow considerations 

There is a considerable body of work in the literature that 
deals with fast dynamo action in simple, steady and three- 
dimensional flows. In our simulation the velocity is chosen 
to have vigorous (fast) dynamo properties in the kinematic 
regime, where the Lorentz force is negligible. One such 
class of flows, which are well known candidates for fast 
dynamo action is the ABC flows. The form of the velocity 
is given by 

uabc = -A(0, sinfcx, cosfcx) +B(cosA;y, 0, sinfcy) (5) 
+ C(sin kz, cos kz, 0) 

This periodic flow is the sum of three steady Beltrami 
waves, parameterized by A, B and C and has the prop- 
erty V x u oc u (velocity is parallel to vorticity). If one or 
more of the constant coefficients (A, B and C) is zero, the 
flow is integrable and is not a fast dynamo. If all three are 
non-zero the flow is non-integrable (Dombre et al. 11986}! 
and contains a mixture of chaotic regions and regular is- 
lands. The special case A = B = C = 1 was introduced by 
Childress H197U}! as a model for kinematic dynamo action. 
Arnold et al. (1983} noted that a steady three-dimensional 
flow with infinite conductivity and chaotic streamlines fa- 
vors the growth of magnetic fields. Their results were fol- 
lowed by the work of Galloway et al. (1984 ) and Moffatt et 
al. (1985) who studied dynamo action in ABC flows with 
finite conductivity. A detailed description of the amplifica- 
tion process responsible for the kinematic dynamo action, 
when the wavenumber of the flow k is equal to one or 
higher, has been provided by Dorch (2000) and Archontis 
et al. H2003[) . The hydrodynamic stability of the ABC flow 
has been studied by Podvigina & Pouquet l}1994fl who 
found that for values of fluid Reynolds number above a 
critical one (Re c = 13) the ABC flow destabilizes, first 
to time dependent but still smooth states, and then to 
a turbulent state. In our simulation, the definition of the 
Reynolds numbers is the same as in Galanti et al. I}1992}l . 

In the following experiments the initial flow is taken 
to be an ABC flow with A = 0.9, B = 1, C = 1.1 and 
wavenumber equal to unity. The reason for choosing val- 
ues of A, B, and C slightly different from unity, is that 
the A = B = C case is too special, because of its exact 
symmetry, and thus turbulence takes longer to develop. 
The purpose of our work is to study dynamo action in the 
turbulent regime of driven ABC flows since astrophysical 
dynamos typically occur in turbulent environments with 
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very high fluid and magnetic Reynolds numbers. However, 
the behavior of turbulent dynamos in the high Re, Re m 
limit may be expected to be to some extent generic. An 
increase in the Reynolds numbers involves smaller and 
smaller length scales where dissipation actually occurs. 
However, the rate of dissipation is determined by the large 
scale developments of the flow and is independent of the 
magnitude of the viscosity once the Reynolds numbers are 
large. Thus, there is good hope that also dynamo proper- 
ties are asymptotically independent of the Reynolds num- 
bers. That hope is indeed borne out by the present numer- 
ical experiments, as is discussed in more details below. 

3. Results 

As a starting point, the forcing function which is imposed 
is such that the ABC flow Eq. © appears. The forcing 
necessary to accomplish this must have the following form: 



f = -^-V 2 u A bc 

tic 



(G) 



However, in order to ensure that the kinetic energy is 
maintained at a value close to the initial value, the am- 
plitude of the forcing is controlled by a differential equa- 
tion, which compensates for the increased dissipation in 
the turbulent phase by increasing the forcing when the ki- 
netic energy tends to drop below the nominal value. The 
amplitude factor in front of the driving force is given by 

-CJkin 

where L is determined from 
dL _ 1 E 

dt Tt.o. Ekin' 



(8) 



with L(t — 0) = 0. The differential equation evolves the 
amplitude factor e L on the turn over time scale, Tt.o. , if 
the kinetic energy tends to deviate from the nominal value 
Eo- The immediate factor Eo/Ekin helps to reduce the time 
delay that results from the integral nature of the exponen- 
tial factor, and the delayed response of the kinetic energy 
on changes in the driving. The kinetic energy control is an 
essential feature, in that otherwise the amplitude of the 
velocity in the turbulent phase would become arbitrarily 
small, at increasingly large fluid Reynolds number, since 
the forcing in the laminar phase needs to scale with 1 /Re. 
It is indeed necessary to maintain the kinetic energy in or- 
der to maintain the same turn over time and actual fluid 
Reynolds number, since the formal fluid Reynolds num- 
ber is defined in terms of the RMS velocity of the laminar 
ABC flow. 

Figure ^ shows the evolution of the kinetic energy in 
time for the first 100 time units. Two phases may be iden- 
tified. The laminar phase which extends up to t — 30 and 
the turbulent regime for t > 30. Note that the kinetic en- 
ergy drops by up to about 40% when the flow enters the 
turbulent phase. After a short adjustment time, the forc- 
ing control returns the kinetic energy to the nominal level, 
where it is then maintained, except for minor fluctuations. 
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Fig. 1. Temporal evolution of the kinetic energy showing the 
transition from a laminar to a turbulent state. 

The magnetic energy for the same time period during 
the laminar phase increases exponentially and the growth 
rate is about 0.15. When the flow enters the turbulent 
phase (£ > 30) the magnetic energy drops slightly and 
later (t — 40) continues to grow exponentially with a sim- 
ilar growth rate as in the laminar phase (Fig. [2J ■ 




40 60 
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Fig. 2. Temporal evolution of the magnetic energy during the 
kinematic regime of the dynamo. 

An analysis of the external work in the above two dis- 
tinct phases shows that the average work in the turbulent 
regime is much larger than in the laminar case (Fig. |3J) . 
On the first hand, the asymptotic level of the average work 
(» 0.25) is independent of Re and incidentally corresponds 
closely to the level that would be obtained for the critical 
value (Re c ). On the other hand, the average work in the 
laminar regime is proportional to 1/Re, and can become 
arbitrarily small for high enough Reynolds numbers. 

This illustrates the point made above, about the forc- 
ing control; without it, the amplitude of the velocity in 
the turbulent phase would drop very much and the kinetic 
energy could not be maintained in order to have constant 
fluid Reynolds number during the simulation. 

The enhancement of the driving force above the lam- 
inar value is shown in Fig. The factor in Eq. (JJJ) is by 
definition = 1 during the laminar phase of the flow but 
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Fig. 3. The evolution of the external average work in time. Fig. 5. The evolution of the viscous dissipation up to t = 100. 



increases during the turbulent phase to maintain the ki- 
netic energy. The asymptotic level of the factor is about 
10 for Re = 100. 
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Fig. 4. The increase of the amplitude of the driving force going 
from the laminar phase to the turbulent phase. 

Figure |S] shows the viscous dissipation as a function 
of time for the laminar and turbulent phase. During the 
laminar phase the dissipation is constant and is balanced 
by the external work. The viscous dissipation increases 
when the flow becomes turbulent. It reaches an asymp- 
totic level when the turbulent velocity is maintained at 
the same level as that of the laminar flow. 

Figure |H1 compares the Joule dissipation to the work 
done against the Lorentz force and also shows the tem- 
poral evolution of the magnetic energy for the kinematic 
regime of the dynamo. 

The work done against the Lorentz force is leading the 
Joule dissipation and is ultimately responsible for the in- 
crease in the magnetic energy. There is a close balance 
between them with a small positive difference during the 
increase of the magnetic energy and small negative dif- 
ference when the magnetic energy decreases, see Fig. Et, 
Fig. Eh- The fluctuations of the Lorentz work lead similar 



fluctuations of the Joule dissipation when the magnetic 
energy increases or decreases. Their difference in the sat- 
urated state of the dynamo fluctuates around zero, as ex- 
pected. The work done by the fluid on the magnetic field 
is, on average, converted into Joule heat and none is left 
to increase the magnetic energy. 

The accuracy in the numerical method is demonstrated 
in Fig. Et, where the dE mag /dt is overplotted (symbols) 
on top of the difference Wl — Qj- These are averages over 
the periodic box, for which the flux term in the magnetic 
energy equation vanish. The magnetic energy equation is 

^p = -V-F poynt -u.f L -Q. T . (9) 

The source of magnetic energy is the work done against 
the Lorentz force and the sink is the Joule dissipation. 
Magnetic energy is transported by the Poynting flux, 
Fpoynt = ExB. According to Eq. ©, dE mag /dt should 
be identically equal to Wl — Qj, and the two are indeed 
very nearly the same. 

Note that the ratio of the Lorentz work to the Joule 
dissipation in Fig. EJl remains close to unity both dur- 
ing intervals of time when the magnetic energy is growing 
and during intervals of time when the magnetic energy is 
decaying. 

4. Magnetic and velocity patterns 

During the laminar phase the flow has points with two- 
dimensional stable manifolds (so called a-type stagna- 
tion points), where the magnetic field is concentrated 
along magnetic "flux cigars" . There are also points with 
two-dimensional unstable manifolds, similar to the /3-type 
stagnation points of the A=B=C=1 flow, where magnetic 
flux sheets are formed (see Dorch|5000). 

Figure[7|(left panel) shows that flux cigars form an in- 
clined triangle around a /3-type stagnation point (located 
at the center of the box). Another one, which is located 
at the upper left corner, is pointing towards the plane of 
the /3-type point while the flow is still laminar. The right 
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Fig. 6. Four panels showing different aspects of the time evolution of the Lorentz work and Joule dissipation during the linear 
regime of the dynamo and comparison with the magnetic energy: a) The Lorentz work Wl (full curve) and the Joule dissipation 
Qj (dashed curve), b) Total magnetic energy E mag . c) The difference between the Lorentz work and the Joule dissipation, d) 
Ratio between work and heating. 



panel of the same figure shows low velocity and high mag- 
netic field isosurfaces at t = 30. The velocity field comes 
into a turbulent phase and the regions of low velocity are 
disturbed, but the triangle formed by the flux cigars is 
still visible. The low velocity isosurfaces no longer contain 
points of exact stagnation, but still correspond to regions 
of vigorous stretching (see Archontis et al. 2003 ) , which is 
why we choose to visualize them. 

Figure is a visualization of the same velocity and 
magnetic field strength for t = 80. The velocity field has 
a more turbulent configuration and the regions with ve- 
locities less that 4% of the peak have disappeared. The 
magnetic flux cigars also disappear and the magnetic en- 
ergy is confined into curved tubes and weaker sheets. 

5. Saturation of the dynamo 

Two experiments were performed to show the dependence 
of the saturation level of the turbulent dynamo on the 
magnetic and fluid Reynolds number. Cases with Rc m = 
Re = 200 and Rc m = Rc = 400 are studied and the 
temporal evolution of the magnetic energy is followed in 
the kinematic and saturated regime. 



Figure shows the growth of the magnetic energy in 
the beginning of the kinematic regime for Re m = Re = 
200 and Re m = Rc = 400. The numerical resolution used 
were 80 3 and 160 3 , respectively. After an initial transient 
phase the growth rate for the second experiment is very 
similar to the growth rate of the first one. The difference 
in the magnetic energy amounts to a factor of 2.7. Thus 
if the magnetic energy for the J60 3 experiment is divided 
with 2.7, it falls right on top of the 80 3 experiment with 
Re m = TOO (see Fig.0. 

Keeping the same values of magnetic and fluid 
Reynolds number, the temporal evolution of the magnetic 
energy is examined in the saturated regime. The mag- 
netic energy eventually stops growing exponentially and 
saturates. A result of considerable significance is that for 
higher values of the fluid Reynolds number, but constant 
magnetic Prandtl number (Pr m = vjrj), the saturation 
level of the turbulent dynamo increases, rather than de- 
creases (Vainshtein & Cattaneo i992). 

Figure shows the temporal evolution of the mag- 
netic energy over a small epoch of the non-linear regime, 
illustrating the effect of increased numerical resolution. 
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Fig. 7. Isosurfaces of magnetic field strength (light) and low velocity (dark) at t = 20 (left panel) and t — 30 (right panel). The 
isosurface level of the magnetic field strength is 75% of the peak value in the snapshot. The velocity isosurfaces correspond to 
4% of the peak velocity in each snapshot. 




Fig. 8. Isosurfaces of magnetic field strength (light) and low 
velocity (dark) at t = 80. 



There is an initial transient, as the larger resolution al- 
lows field concentrations to collapse further. Note that 
there arc then sections of time where the growth rates 
and the decay rates are quite similar in the Re m = 400 
and Re m = 200 cases. The sections after the largest maxi- 
mum in the two cases, and the rise towards the next max- 
imum are two examples. This might be fortuitous, but 
more likely it corresponds to parallel evolution of similar 
magnetic structures. 
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Fig. 9. The growth of the magnetic energy in the kinematic 
regime when Re m = Re = 200 (straight line) at 80 3 resolution 
and Re m = Re = 400 (dashed line) at 160 3 resolution. The 
latter divided by 2.7 is shown by the lower (dashed) line. 

The overall level is slightly higher, and the evolution 
is slightly delayed in the high Re case, because flux struc- 
tures that were barely resolved at 80 3 are allowed to col- 
lapse in the 160 3 case. This is indeed also what happens 
in the kinematic case, but here the size of most structures 
is controlled by the Lorentz force, and not by diffusion, 
so the effect is only marginal. For the same reason, the 
effect may be expected to vanish altogether as the mag- 
netic Reynolds number goes to infinity. In addition, the 
size of the magnetic structures in the kinematic regime 

1/9 

scales as Rc m ' . In the saturation regime the size is again 
determined by the balance between the advection and dis- 
sipation of the magnetic energy and presumably changes 
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following the same scaling with no crucial dependence on 
the value of Re. 
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Fig. 10. Temporal evolution of the magnetic energy in the 
saturated regime for Re m = 200, Re = 200 (solid line) and 
Re m = 400, Re = 400 (dashed line). 



After some time, the saturation level for Re m = 400, 
indeed fluctuates around a level that is only marginally 
larger than for Re m = 200. This is small compared to the 
factor 2.7 for the kinematic regime of the dynamo. 

One also notices that there are more fluctuations in 
the higher resolution case (e.g. the extra bump on the way 
down after the largest maximum), as a result of the small 
scales resolved when going to higher Re m . That trend is 
expected to continue with increasing magnetic Reynolds 
number. 
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Fig. 11. Temporal evolution of the magnetic energy for Pr m > 
1. Numerical resolution is 96 3 for Re m = 300 and 80 3 for Re m = 
200. 
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6. Magnetic Prandtl number and growth rate 

In Section [3] we examined the growth rate and the satu- 
ration level of a fast dynamo when the magnetic Prandtl 
number is equal to one. Equally important is the study of 
the degree of their dependence when Pr m is varied. Thus, 
additional experiments were performed where the forcing 
again is such that it maintains the kinetic energy at a value 
close to the initial one, apart from small oscillations. 

Not surprisingly, the growth rates during the lami- 
nar phase (t < 30) of the flow are found to be the same 
(Figures If II I12JI . By construction, the velocity is exactly 
the same in the laminar phase, independent of Re, and the 
effect of increasing Re m is simply, as in the comparison of 
the Re,„ = 200 and Re m = 400 experiments, to give a 
different transient at the beginning of the simulation. 

Less trivially, the average growth rates in the turbu- 
lent regime also turn out to be similar with no crucial 
dependence on the magnetic Prandtl number for high Re. 
This is illustrated in Figure where the difference in 
the magnetic energy during the turbulent phase (t > 30) 
is very small between the two experiments with Re = 200 
and Re = 300. This may also be seen as a consequence 



Fig. 12. Temporal evolution of the magnetic energy for Re m = 
200 and Pr m < 1. 



of the generic properties of high Reynolds number tur- 
bulence. More specifically, by one of Kolmogorov's main 
assumptions, energy dissipation (and hence the forcing) 
is independent of the fluid viscosity once the Reynolds 
number is high enough. Thus, and by analogy with the 
previous experiments, we expect that the only important 
change in the magnetic field is the development of more 
fine structure with increasing magnetic Reynolds number. 
The growth rate is determined by the large scale proper- 
ties of the flow, and in particular by the turn over time, 
which by construction does not change between these ex- 
periments. 

Finally, the above experiments reinforce the conclusion 
from Section |SJ that the saturation level of the dynamo 
does not depend crucially on the value of Pr m , but only 
increases slightly with Re m as it is shown in Figure ITT1 
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7. Nature of the dynamo 

The results of the above non-linear experiments are sim- 
ilar to the results of previous studies of turbulent hydro- 
magnetic convection (Nordlund et al. I1992|l and both re- 
sults indicate that during the dynamo action there is a 
close balance between the work done against the Lorentz 
force and Joule dissipation with a relatively small positive 
(negative) difference when the magnetic energy increases 
(decreases). In other words, work and dissipation are in 
near balance, with a relatively small difference responsi- 
ble for the average growth of the magnetic energy. One 
may rightly ask why this is so; is the small positive dif- 
ference then just fortuitous, and might the difference not 
just as well be negative for a slightly different flow? 

To address this question we make use of a kinematic 
dynamo experiment, where the analysis is simpler than 
in a turbulent case. Turbulence makes the analysis of 
the Lorentz work difficult — there is much time dependent 
noise, which corresponds to waves averaging out over time. 
The work is much less noisy when the velocity is prescribed 
and/or has not yet entered the turbulent phase, and we 
therefore analyze the results of an experiment that deals 
with kinematic dynamo action. The velocity field is cho- 
sen to be the A = B = C = 1 flow with k = 1 and the 
magnetic Reynolds number is Re m = 300. 

The convergence and the diffusion rates at the points 
of strong flux concentration are of the order of y2, which 
is much larger than the dynamo growth rate (which is 
about 20 times smaller, e.g. see Childress & Gilbert[1995), 
so the advection and diffusion are in close balance in the 
neighborhood of the "flux cigars" . During the exponential 
growth of the magnetic energy the size of the flux cigars 
does not change dramatically, but the magnetic energy 
increases gradually. 

A detailed analysis for the latter case confirms that this 
is indeed the case, and that the bulk part of the Joule dis- 
sipation and the work done against the Lorentz force are 
in detailed balance. This is understandable as a balance 
between kinetic energy being converted into magnetic en- 
ergy, that is then immediately converted into heat through 
Joule dissipation. 

Significant net positive work (defined by the difference 
Wl — Qj) occurs instead in regions where the magnetic 
field is weak, and the dissipation is very weak as well. In 
these regions, which cover 90% of the volume in this ex- 
ample, there is almost "pure work" ; the Joule dissipation 
is much smaller than the work, and the net work is almost 
identical to the net work in the whole box. 

In other words: The dissipation is strong only in the 
neighborhood of the ev-type stagnation points, but there 
it is balanced almost exactly by advection/convergence. 

Table shows the estimated values of the average 
work and dissipation over the 89% of the computational 
volume. On the one hand, the dissipation is very weak (less 
than 0.12% of the maximum dissipation level) and on the 
other hand, almost all the net work occurs in those 89% of 
the volume. Moreover, the largest amount of dissipation 



Average Lorentz work 


326.3 


Average Joule dissipation 


271.1 


Average magnetic energy 


332.3 


Energy growth rate 


0.166 


Fractional dissipation level 


0.0012 


Fractional volume 


0.89 


Net work in weak dissipation regime 


55.1 


Net work in high dissipation regime 


0.09 



Table 1. Summary of average values from the kinematic ex- 
periment with Re m = 300. 



and work occur in the remaining 11% of the volume and 
they are in almost perfect detailed balance. 

The above results lead to the conclusion that dynamo 
action occurs primarily in regions where the field is weak. 
The net work is much higher than the dissipation and as a 
result the field is amplified by dynamo action. Most of the 
net work occurs when the weak field is bent and stretched. 
The magnetic energy is then distributed over the volume 
through the Poynting flux and the field lines are dragged 
out and pile up against the local flux cigars where balance 
occurs between stretching and diffusion. The weak field is 
almost perfectly advected with the fluid and this advection 
of the weak field would be possible even in the saturated 
regime where the flow still has a good grip on the weakest 
field, although the strongest field is now exerting a strong 
resistance against motion/stretching, through the Lorentz 
force. 

The above picture of the dynamo process is most likely 
related to the observation that the growth rates found for 
the kinematic ABC dynamos are always consistent with 
the large scale turn over times, even though the ratio 
of the magnetic energy to the rate of Joule dissipation 
(Emag/Qj), which is the time over which the magnetic 
energy is destroyed and replenished, is much shorter. 

In conclusion, the work done on the weak part of the 
field is responsible for the increase in the magnetic energy 
by bending and stretching the magnetic field lines. The 
exponential amplification dominates the Joule dissipation 
which is weak through almost all of the volume. The de- 
tailed balance between the Lorentz work and the dissipa- 
tion occurs only around the points where strong magnetic 
field structures are formed. Thus, a dynamo appears to be 
a process where amplification of the weak field plays the 
dominant role in the growth of the magnetic energy. 

8. Conclusions 

In this paper turbulent dynamo action has been studied 
and some new aspects of kinematic dynamo action have 
been illustrated. It was found to be important to keep the 
kinetic energy level of the same order through the sim- 
ulation, by automatically adjusting the amplitude of the 
driving. Using relatively high values of Re m and Re we 
find that the growth rate of the magnetic energy amplifi- 
cation apparently depends on the stretching ability of the 
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flow rather than on the values of the above parameters. 
This is obvious even in the turbulent phase, for different 
values of the magnetic Prandtl number, where the dissipa- 
tion is increased. The level at which the magnetic energy 
saturates increases slightly with Re m but the transient in 
the saturated regime is expected to become insignificant 
at sufficiently high Re m . 

An important result is the existence of a near bal- 
ance between the Lorentz work and the Joule dissipation, 
which is apparent in both the laminar and the turbulent 
regime of the forced ABC flow. Remarkably, the balance 
originates primarily from small regions where strong mag- 
netic flux structures are concentrated. The nature of the 
dynamo is such that the net growth of magnetic energy 
comes about through stretching and folding of the weak 
magnetic field in the rest of the volume, where dissipation 
is very weak. 

The aim of this paper was not to study dynamo action 
on a global scale (including aspects such as differential 
rotation, stratification and convection) but rather to study 
the different stages of dynamos produced by a family of 
hydrodynamically unstable flows. However, the physical 
processes and the dynamo mechanisms studied here are 
likely to be of similar importance in more realistic stellar 
dynamo settings. 
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